Topological superconductivity in bilayer Rashba system 
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We theoretically study a possible topological superconductivity in the interacting two layers of Rashba sys- 
tems, which can be fabricated by the hetero-structures of semiconductors and oxides. The hybridization, which 
induces the gap in the single particle dispersion, and the electron-electron interaction between the two layers 
leads to the novel phase diagram of the superconductivity. It is found that the topological superconductivity 
without breaking time-reversal symmetry is realized when (i) the Fermi energy is within the hybridization gap, 
and (ii) the interlayer interaction is repulsive, both of which can be satisfied in realistic systems. Edge channels 
are studied in a tight-binding model numerically, and the several predictions on experiments are also given. 
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The topological aspects of the electronic states in solids 
have recently attracted intensive interest. In addition to the 
quantum Hall effect |fl] [2), the anomalous Hall effect Q, 
spin Hall effect [4], and topological insulators (TI's) J5] |6) 
turn out to be topological phenomena driven by the Berry 
curvatures JT| in momentum and/or real spaces. Especially, 
the topological insulators are realized in the time-reversal 
(T) symmetric systems with the relativistic spin-orbit inter- 
action (SOI), and are characterized by the Z2 topological in- 
dices. On the other hand, the concept of the topological super- 
conductors (TSC) or superfluidity (TSF) has been proposed in 
the context of 3 He 0. Recently the p-wave superconductiv- 
ity with broken T symmetry in Sr 2 Ru04 has been discussed 
from this respect ||9] [lOl . Helical topological superconduc- 
tors with T symmetry have also been discussed ATI . On the 
other hand, the unified scheme of the classification of the TIs 
and TSCs according to the three symmetries (T symmetry, 
particle-hole symmetry, and chiral symmetry) and the dimen- 
sion of the system has been established fi"2l . 

The implementation or fabrication of the TSCs is an im- 
portant issue especially because it is expected to support the 
Majorana fermions lfT3l - IT7l . a promising tool for the quan- 
tum information processes. An interesting recent proposal 
is the proximity-induced superconductivity of the surface 
state in the three dimensional TI lfT8l . Because of the elec- 
tron fractionalization in the surface Dirac fermions, Majorana 
fermions are expected to appear at the interface between the 
ferromagnet and s-wave superconductors both of which are 
on top of TI. It is also proposed that the superconductivity in 
a doped TI Cu r Bi2Se3 lfT9ll can be topological [ 20-24 1. Be- 
sides, the possible TSC in noncentrosymmetric superconduc- 
tors with the Rashba spin splitting has been discussed ||25ll26l . 
The basic idea is that the spin direction rotates with the mo- 
mentum for the spin split bands, leading to the p ± ip pair- 
ing form when the superconducting order parameter is pro- 
jected onto each of the bands. However, it has been clari- 
fied that due to the interference between the two order pa- 
rameters in the two spin split bands, spin triplet p-wave pair- 



ing should be dominant over spin singlet s-wave pairing for 
the realization of TSC l|25l l27l |28l . One possible way to re- 
solve this difficulty is to delete one of the spin split bands 
by external magnetic field or by the exchange splitting due to 
the proximity to a ferromagnet, leaving alone the p + ip (or 
p — ip) pairing with broken T symmetry ll26l [29 - 33 1 anal- 
ogous to the case of Sr 2 Ru04. Experimentally, it is known 
that the superconductivity emerges at the interface of LaA103 
and SrTi03 |34|, and the role of Rashba interaction has been 
identified there [35]. Therefore, the two-dimensional super- 
conductors with Rashba coupling are now available. It has 
been also reported that the superlattice of a heavy-fermion 
compound, CeCoIn 5 , can be fabricated and its superconduc- 
tivity is observed 11361 . Theoretically, a robust Fulde-Ferrell- 
Larkin-Ovchinnikov state has been proposed in the presence 
of SOI [37 1 ■ Also, up to now, there have been theoretical pro- 
posals about interface induced superconductivity [38 39] and 
superconductivity in multi layer system |40|. The advantage 
of the interface or heterostructure systems is the great con- 
trollability of the structure, symmetry, doping, and SOI. For 
example, the bilayer quantum well to form the two interacting 
two Rashba system can be fabricated both in the oxides fiTll 
and semiconductors |42|. 

In this Letter, we theoretically study the superconductivity 
in the interacting two layers of Rashba systems. In contrast 
to the previous proposals to delete one of the spin split bands, 
this system doubles the degrees of freedom without breaking 
T symmetry, giving richer possibilities. The model system 
we consider is schematically shown in Fig.[TJa). At each inter- 
face, the two-dimensional electron gas (2DEG) is formed, and 
this sandwich structure produces two 2DEG's with the tunable 
interaction by the width of the middle layer. There are two 
types of the interaction. One is the transfer or hybridization 
of the electrons, and the other is the electron-electron interac- 
tion. The former induces a gap in the band structure, denoted 
by e. The Hamiltonian for the kinetic energy is given by 

k 2 

Hq = —I s ® I<r - £l s ®o- x + a (k x s y - k y s x ) a z , (1) 
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FIG. 1. (color online), (a) (left panel) Schematic view of the model 
with the interfaces between two kinds of materials A and B. An ex- 
ample is that A is LaA103 and B is SrTiOs. 2DEGs are formed in 
the SrTiOs side and we expect these two hybridize with sufficiently 
short width of the middle layer, (right panel) Schematic wave func- 
tions of bonding and antibonding states along the z direction. The 
hybridization splits two Rashba bands into bonding and antibonding 
states with the hybridization gap e. (b) Dispersion of the Hamiltonian 
of bilayer Rashba model Eq. QJ. The position of the Fermi energy 
(A-C) is crucial to nontrivial phases. The phase diagrams for each 
case is shown in Fig. [2] 



where s and a respectively denote the Pauli matrices for spin 
and layers, and a represents the strength of Rashba SOI. The 
eigenvalues of this Hamiltonian is 



E k = k 2 /2m ± Ve 



2 k 2 = k 2 /2m±X k . 



(2) 



Note that these two layers are exposed to the local electric 
field with opposite directions, and this is described by a z in 
the 3rd term of Eq. The bands shown in Fig. [TJb) are 
doubly degenerate at each fc-point in the whole Brillouin zone 
because this system has both T symmetry, T — is y K, where 
K is the complex conjugation operator, and inversion (I) sym- 
metry, I — <j x . We assume the following density-density in- 
teraction between the electrons 

(x) = -U (nj (x) + n\ (a?)) - 2Vn x {x) n 2 (x) , (3) 



where n n 



=i,2 = 4. c i7s c <?s is the electron density in 

layer 1 and 2. U and V are intralayer interaction and in- 
terlayer one, respectively. We construct the Bogoliubov de 
Gennes (BdG) Hamiltonian with mean field approximation 



H 



BdG : 



dfc* 



"Ho (k) - n 
A 



A 

-(Ho (k) - /i) 



(4) 



r is the Pauli matrices in Nambu (particle-hole) 
space and the basis are in the form of ^J, = 

( c ifcT' c 2kV c i*4' c 2*4' C i-H) C 2-H, -ci-fef, -c 2 -kt)- 
Note that the BdG Hamiltonian Eq. Q and the analysis 
below is similar to that in Ref. EU1 . However, the situation 
we consider is quite different from theirs, i.e., (i) we consider 
the 2D case while Ref. GUI considered the 3D case, (ii) the 
energy dispersion in Eq. (EJ of the Rashba system is different 
from that of Dirac fermion due to the presence of the k 2 terms 
in the kinetic energy, (iii) there is a tunable parameter a in our 
model, and (iv) we consider the layer-index while Ref. l20l 
considered the orbital index, and the meaning of U and V is 
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TABLE I. There are four possible nonvanishing pairing potentials in 
our model with the assumption Eq. {3}. Ai, A2, A3, and A4 be- 
long to A\ % , Aiu, A2U, and E u irreducible representations of D^h, 
respectively 1 43 1 . Matrix representations are off-diagonal elements 
of BdG Hamiltonian, i.e., entries in the column "explicit representa- 
tions" are products of the corresponding entries in the column 
"matrix", r x and ^. 



different. We consider pairing potential A (k) according to 
the lattice structure and the interaction Hamiltonian Eq. 0). 
Since we are conscious of the interface between SrTi03 and 
LaA103, the lattice symmetry is assumed to be D^h , which 
contains following operations: a four-fold rotation about z 
axis, a twofold rotation about x axis and a mirror reflection 
with respect to xy plane. Note that the mirror reflection sends 
a site on a layer to that on the other. In the weak coupling 
limit with purely short-range interaction, k independent 
pairings are favored compared with k dependent anisotropic 
pairings. Only 6 forms listed in Table [I] can have nonzero 
values among the 16 possible products of (1, s x , s y , s z ) 
and (1, & X) <j yt (j z ). Ai is the even combination of the 
intralayer pairings {c\^ci^-\- 024-024.) and the interlayer pairing 
(C14-C24 — 014024-) , and is parity even. It is remarked that the 
parity operation here represents (r,a) — > (— r, o) with the 
layer index a (a represents the opposite layer to a), and is not 
the even-odd symmetry with respect to the exchange of the 
two coordinates (n, r 2 ) —> (t*2, 1*1) of the pairing electrons. 
A2 is the interlayer pairing with odd parity, while A3 is the 
intralayer with odd parity. Last, A4 is the interlayer pairing 
with odd-parity. 

The excitation energy of quasiparticles are obtained by di- 
agonalizing the BdG Hamiltonian Eq. (4) with fixing the pair- 
ing potential to each Aj (see note |44|). We find supercon- 
ducting gap for A4 has point nodes (in the k x direction when 
one chooses s y a y ), and the others have full gap. We estimate 
the superconducting critical temperature Tq by analyzing su- 
perconducting susceptibility for each pairing potentials. The 
pairing susceptibility \o is defined as 



Xo 



-T^^Tr^Go (k)T x G (fe)] 



l-2/(& + A fc ) , l-2/(&-A fc ) 



2(& + A fc ) 



2 (Cfc — Afe) 



,(5) 



where f(E) is the Fermi distribution function and = 
k 2 /2m — fi. Hereafter, we assume the weak coupling limit. 
The other susceptibilities can be calculated by replacing t x 
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with t x s<j (sa is a matrix representation for each pairing in 
Table[l]). As a result, they can be expressed by xo> which con- 
tains the logarithmic divergence estimated by the value at the 
Fermi surface. The linearized gap equations are obtained as 
(two of them are shown explicitly) 



case: A 



U\o - 1 Uxoi 
Vxio Vxi - 1 



u + v{e/x k y 



V X 2 = V 



i - (e/x k y 



Xo-l=0 for Ax, (6) 
Xo = l forA 2 , (7) 



when the Fermi surface lies in the hybridization gap (case 
B). Note that, for Ai, there are two choices of vertex, t x I 
or r x a x . Correspondingly, the gap equation has 2x2 form 
with the coefficients xo — ~ / d£D + (ri)(t&nh((3ri/2) /2r/), 
Xoi = Xio = (e/A fc )xo and xi = (e/X k ) 2 Xo- In the lin- 
earized gap equations for A 3 and A4, the susceptibilities are 
X3 = 2x4 = [1 — (£/Afe) 2 ]xo an d the coefficients are U and 
V, respectively. \ k is evaluated at the momentum crossing 
the Fermi energy since we assume the weak coupling limit. 
D + {r]) is the density of states of the outer cone in the band 
structure [Fig.[TJb)] and 77 is energy measured from the Fermi 
energy. From the highest T c , we obtain the phase diagram in 
the UV plane as shown in Fig. [2] (case B). It is seen that the 
phase diagram consists of Ai, A2, A3, and nonsuperconduct- 
ing regions. A4 cannot be the leading instability in this phase 
diagram, and hence all the pairing states are fully gapped. 
The phase boundaries depend only on U/V since they are de- 
termined by the comparison among the effective interactions 
for each pairing, which contain terms linear in U and/or V. 
When both U and V are repulsive, the superconducting in- 
stability is absent, while some pairing occurs in all the other 
cases. When the intralayer attraction U{> 0) is dominant, the 
conventional pairing Ax occurs, while the unconventional A 2 
and A3 are realized in other regions. Starting from the single 
layer case where the superconductivity is realized as in the 
case of LaA10 3 /SrTi03 interface 11341 . U is considered to be 
attractive. Assuming that this attractive interaction U is given 
by the short-range electron-phonon coupling reduced by the 
Coulomb repulsive interaction, the repulsive interlayer |V| is 
expected to be larger than U. Therefore, we expect the situa- 
tion where A3 is realized in the bilayer Rashba system. Note 
that although the region of A3 is independent of the Rashba 
SOI, the effective coupling constant for A 3 is U[l — (e/X k ) 2 ] 
which is zero for a = 0. Therefore, A3 is the superconduc- 
tivity induced by Rashba SOI [45]. When there is no SOI, 
the bands are simply bonding and antibonding ones with the 
hybridization gap, where the interactions U and V drive only 
conventional spin singlet pairing. The spin states are mixed 
with finite SOI, and this causes unconventional spin triplet 
pairing. According to the symmetry of pairing, U drives Ax 
and A3 while V drives Ax, A 2 and A4. These features are 
valid for both cases with single and multi Fermi surfaces. Ad- 
ditional Fermi surface in case A and C gives another logarith- 
mically divergent term and modifies phase boundaries. In case 
A, A3 instability is always weaker than that of Ax, and A 2 in- 
stability is also weaker when SOI is less than a certain critical 




FIG. 2. (color online). Phase diagram of superconductivity in the 
UV plane. Yellow, brown, and orange indicate the region where Ax, 
A2, and A3 show the strongest instability, respectively. Here we 
take the unit where m is the unity. The white region is nonsupercon- 
ducting. Each row corresponds to the Fermi energy at A, B, and C in 
Fig.[TJb). The system can be topological only when the Fermi energy 
lies in the hybridization gap, i.e., case B, and in this case A2 region 
expands as the SOI increases while A3 region is independent of the 
strength of SOI. However, SOI is indispensable for A3. See the text. 



value. In case C, the boundary between Ax and A3 is deter- 
mined by the strength of SOI. Note that case C can appear 
only when SOI is strong enough for the band structure to have 
double minima [see Fig. [TJt>)]. There can be rapid change 
of phase boundaries between case A(C) and B because the 
density of states in 2D has step-function-like behavior at the 
edge of band bottom or top, which gives an additional log- 
arithmic singularity. Apart from topological nature, the fact 
that various kinds of pairing potential occur corresponding to 
the strength of SOI is an interesting result in this Letter. 

Next we consider the topological nature of each pairing 
state. According to the generic classification scheme fl2l . the 
present system belongs to Dili in 2D. Therefore, the homo- 
topy class is characterized by Z 2 number in contrast to the 
Z in 3D. We can classify the superconductors into two cate- 
gories, i.e., topologically trivial superconductors and the he- 
lical one. Since our system has T symmetry, particle-hole 
symmetry, and I symmetry in the normal state Hamiltonian 
Eq. ([TJ, the criterion in Refs. 112011461 can be applied. Namely, 
the present system becomes a topological superconductor if 
(i) it has odd-parity pairing symmetry with full superconduct- 
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ing gap (A2 or A3) and (ii) odd number of Kramers pairs at 
the time-reversal invariant momenta T Q 's are below the Fermi 
surface. The condition is explicitly written in 



(-ir=n(-!) 



JV(r„) 



(8) 



case: A 



case: B 



1H1 







case: C 
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where N(T a ) is the number of occupied Kramers pairs at T a . 
To satisfy these conditions in our system, the Fermi surface 
must lie in the hybridization gap, i.e., case B. In terms of the 
symmetry of the pairing potential, there are two possibilities, 
A 2 or A 3 , for TSC. 

Next, we confirm the system is exactly in the topologi- 
cal phase under such conditions. It is known that a nontriv- 
ial topological number in a bulk state accompanies topologi- 
cally protected states localized at an edge in a finite system - 
bulk-edge correspondence. We study following tight-binding 
model describing continuous model Eq. Q in the low energy 
regime 

H = ^d T t rT ,c r i+ l^4A4-h.c. I -^(// + ea x )clc r , 

rr' \ r / r 

(9) 



\—t±it's v a z for r — r'±a x e x 
t rr , = { v . (10) 

I— t^fivs x a z for r = r' ±a y e y 

The hopping parameters are related to the ordinary param- 
eters in Rashba system in the form of t = — l/2m, t' = a/2 
and \j! = fi — 2 jm. We consider the A3 pairing state in the 
cylindrical sample, i.e., open boundary condition to the x di- 
rection and periodic one to the y direction. The energy spec- 
trum of this model is shown in Fig.[3]with different positions 
of the Fermi surface. One can see that there are states cross- 
ing bulk superconducting gap only when the Fermi surface 
lies within the hybridization gap (case B). The right-going 
and left-going states are doubly degenerate, respectively, and 
localized at opposite edges. These are nothing but the heli- 
cal Majorana edge modes. Then, we can conclude that when 
the bulk has nontrivial Z 2 number, corresponding helical edge 
states appear at the edge of the system. Helical edge states are 
also generated in the case of A2. However, there are no edge 
states with Ai. These results are totally consistent with the 
topological nature of the bulk states. 

This helical superconductivity has several unique features 
compared with the chiral one [25 47 H491 . Since the T sym- 
metry is preserved, there remains the Kramer's degeneracy for 
the states. Therefore, the spin degrees of freedom is still alive 
for the Majorana edge channels, and Majorana bound state 
in the core of the vortex. This fact leads to the several non- 
trivial physical consequences. One possible experiment is the 
Andreev reflection at the interface of the normal metal and the 
helical superconductor [25 1. Because of the two counter prop- 
agating Majorana edge channels with opposite spins, the spin 
conductance is strongly dependent on the incident angle of 



FIG. 3. (color online). Energy spectrum of the cylindrical sample 
with A3 pairing for different positions of the Fermi surface, i.e., case 
A, B and C (from left to right); see Fig.fJJb). These figures show 
that helical edge modes appear at each edge of the sample in the 
superconducting gap only when the Fermi surface lies within the hy- 
bridization gap. 



the electrons, and also very sensitive to the external magnetic 
field B. Since B acts mostly through the Doppler shift, even 
a tiny B can induce the large spin conductance of the order of 
the charge conductance ESI . Magnetic impurity at the helical 
edge channels is another interesting possibility [47|. Because 
of the Majorana nature, the spin components are reduced from 
3 to 1, i.e., it is of Ising-like. This fact leads to the novel 
Kondo effect of the magnetic impurity at the edge of a helical 
superconductor showing the anisotropic and strong temper- 
ature dependent magnetic susceptibility. Josephson junction 
between the two helical superconductors offers a laboratory 
for the even more rich physics 11481 . Because the forward 
scattering channels are open with the 4 species of Majorana 
channels, Tomonaga-Luttinger liquid can be realized in the 
two interacting sets of the helical edge channels at the Joseph- 
son junction. This leads to the novel temperature and bias de- 
pendence of the quasiparticle tunneling conductance through 
this junction l48l . Last, the appearance of the two zero en- 
ergy Majorana bound states at the core of the vortex is ex- 
pected i50ll5T1 . but the detailed analysis of this subject is left 
for future studies. 

In summary, we have studied the possible topological su- 
perconductivity in the interacting bilayer Rashba systems, 
which can be realized at the interfaces of oxides and semicon- 
ductors. When the Fermi energy is in the hybridization gap 
and the interlayer interaction is repulsive and stronger than 
the attractive intralayer interaction, the topological supercon- 
ductivity is realized. The helical Majorana edge channels are 
confirmed numerically in a tight-binding model in this case. 
Andreev reflection, Josephson junctions, and the Kondo ef- 
fect are the possible phenomena to demonstrate the unique 
features of the helical Majorana modes which have spin de- 
grees of freedom. 
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